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Abstract 

Nonlinear dynamics of the thermal and electromagnetic instabilities of the mixed state in type II 
superconductors has been analysed taking into account the effect of dissipation and dispersion. The 
existance of nonlinear running waves describing the final stage of evolution of the thermomagnetic 
instability in superconductors is demonstrated analitically. 

The evolution of the thermal (T) and electromagnetic {E^ H) perturbations is described 
by the nonlinear equation of the thermal conductivity [1] 

dT 

iy^ = V[^VT]+jE, (1) 

by the Maxwell equations 

IdH 

rotE = (2) 

An ^ 

rotH = —J (3) 

and by the equation of the resistive state 

j=l{T,H)+liE), (4) 

where u = z^(T) and k = k{T) are the heat capacity and the thermal conductivity respec- 
tively; jc is the critical current density and jr is the resistance current density. 

The above system is essentially nonlinear because the right-hand part of Eq.(l) contains a 
term describing the Joule heat evolution in the region of the resistive phase. Such a set (1)- 
(4) of nonlinear parabolic differential equations in partial derivatives has no exact analytical 
solution. 

Let us consider a planar semi-infinite sample (x > 0) placed in external magnetic field 
-> dH 

H = (0, 0, He) growing at a constant rate = const. According to the Maxwell equation 

(2), there is a vortex electric field E = {0,Ee,0) in the sample, directed parallel to the 
current density j: E \\ j; where H^, is the amplitude of the external magnetic field and E^ is 
the amplitude of the external electric field. 

The solution of the system of equations (l)-(4) may be presented as a function of new 
automodel variable ^{x,t): 



1 



T^T[C{x,t)], 



(5) 



Substituting (5) into the system (l)-(4) gives, as a result of simple differentation, the 
following system 



dt 



dT 



(P^dT fdjV d'^T 



+ \jc{T)+3r{E)]E, 



(6) 



d'^idE fd^y d^E _4:n 
dx^ d^ \dx I d^"^ 



djc dT djr dE 
dT~d^^dE~d^ 



dt 



(7) 



In order the system (6), (7) was only function from ^ at the substation (5) it is required 
carrying out the following conditions: 



(8) 



(9) 



(10) 



where A, C, G are functions from ^, the type of which will be determined below. Solving 
first two system of equations (8)- (10) we have a relation 



Whence just filows relationship between G and A 

til 



(11) 



(12) 



where m is a free constant of integrating of the equation (11). From (8) and (9) follows that 
^{x,t) must satisfy single-line equation in private derivation 



d^ _ d^ 
dt dx 

the only solution of which is the function 

^{x,t) ^F{x-ut) 
Using (14) we can immediately obtain 



(13) 



(14) 
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= A{0^-Fu, C(0 = 0. (15) 

It is possible to ensure F = 1 by the transformation of coordinates and time. Thereby, we 
find final automodel substitution 

^ — X — ut, (16) 
corresponding to solution of running type wave [2] . 

For the automodelling solution of the form (16), describing a running wave moving at a 
constant velocity v along the x axis, the system of equations (l)-(4) takes the following form 

dT 

-.[A^(T)-7V(To)] = .--— (17) 
dE Anv 



di c 



j, (1^ 



E = -H. (19) 
c 



The thermal and electrodynamic boundary conditinos for equations (17)-(19) are as follows: 

r(e ^ +oo) = To, ^(e ^ -oo) = 0, 

^ (20) 



E{C +oo) = 0, E{C -oo) = E, 



where Tq is the temperature of the cooling medium. 

Let us consider the Bean- London model of the critical state for the dependence jciT, H) 
[3] 

je(r)=jo[l-a(T-ro)] (21) 

where jo is the equilibtium current density, a is the thermal haet softening coefficient of the 
magnetic flux pinning force. The characteristical field dependence of jr{E) in the region of 
sufficiently strong electric field {E > Ef) can be aproximated by the piecewise linear function 

2 

jr ~ CTfE, where (7/ = ~ GnHc^/H is the effictive conductivity in the fiux flow regime; 

7] is the viscous coefficient, $o = is the magnetic flux quantum, (T„ is the conductivity in 

the normal state; Ef is the bouncFary of the linear area in the voltage-current characteristics 
of the sample. 

Excluding variables T(^) and H{^) from Eqs.(17) and (19) , and taking into account the 
boundary conditions (20), we obtain an equation describing the electric fleld E{^) distribu- 
tion (£^-wave): 

d^E 



de 



+ 



Aiwdj^dE AiTv^a N{T) - N{To] 
~^dE~di ^ ~^ ;<T) 
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where the dependency T = T (e, defined by expression (2), (4) and have the form 



T 



T.,'- 

a 



, A.'KV dE 



(23) 



Here N{T) ^ Jv{T)dT. 



dE. 



The analysis of the phase plane {E, --^) of Eq.(22) shows that there are two equilibrium 

points: Eq = 0, T = Tq is the stabile node and E = Eg, T = T* = T{Ee,0) is the saddle. 
The solution may be rcprcscntedof in the form of the shock- wave-type with amplitude Ee, 
by joining these two equilibrium points. The velocity of E-wsive is determined by the Eq. 
(22) with account of the boundary conditions (20): 

c2 E^ 



Stt 



N 



1 



To + -\jc{T)+jrm 



N{To) 



(24) 



Using Eqs. (19) and (22) we find the expression for distribution of the magnetic field H in 
the case of iJ-wave : 



'de^~^ 



djr 



dE 



E= 



dH N{T)-N{To) 

c 



av H"^ 



0. 



(25) 



The velocity of H wave is connected with its amplitude by the following expression 

tt2 

N{T) - N{To) = (26) 



The condition (26) reflects the adiabatical character of the wave propagation: the magnetic 
fleld energy transferred by the wave ensures local heating of the sample in the close vicinity 
of the wave front. Thereby, depending on the external conditions at the sample surfase there 
may exist two kinds of thermomagnetic waves. 

The system of Eqs. (l)-(4) is invariant with respect to an arbitrary translation. There- 
fore, the wave propagation conditions can be found for an arbitrary critical current density 
depending on T and H. The results can also be obtained for an arbitrary temperature 
dependence of thermophysical parameters ly and k of superconducting material and for an 
arbitrary function jr{E). 



Reference 



[1] N.A.Taylanov Superconduct. Science and Technology, 14, 326 (2001). 



4 



[2] Karpman V.I., Nonlinear Waves in Dispersive Media (Nauka, Moscow, 1973). 
[3] Bean CP. Phys. Rev. Lett. 8, 6, 250(1962). 



5 



NIZAM A.TAYLANOV 



Address: 

Theoretical Physics Department, 
Institute of Apphed Physics, 
National University of Uzbekistan, 
Vuzgorodok, 700174, Tashkent, Uzbekistan 
Telephone: (9-98712)-461-573 
fax: (9-9871)- 144-77-28 
e-mail: taylanov@iaph.tkt.uz 



6 



